Trapped one-dimensional ideal Fermi gas with a single impurity 
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Properties of a single impurity in a one-dimensional Fermi gas are investigated in homogeneous 
and trapped geometries. In a homogeneous system we use McGuire's expression [J. B. McGuire, 
J. Math. Phys. 6, 432 (1965)] to obtain interaction and kinetic energies, as well as the local 
pair correlation function. The energy of a trapped system is obtained (i) by generalizing McGuire 
expression (ii) within local density approximation (iii) using perturbative approach in the case of a 
weakly interacting impurity and (iv) diffusion Monte Carlo method. We demonstrate that a closed 
formula based on the exact solution of the homogeneous case provides a precise estimation for the 
energy of a trapped system for arbitrary coupling constant of the impurity even for a small number 
of fermions. We analyze energy contributions from kinetic, interaction and potential components, 
as well as spatial properties such as the system size. Finally, we calculate the frequency of the 
breathing mode. Our analysis is directly connected and applicable to the recent experiments in 
microtraps. 

PACS numbers: 03.75.-b,67.85.-d 
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Ultracold Fermi gases in the uniform space and in 
harmonic traps have been studied both theoretically[l] 
and experimentally [2], exhibiting fascinating phenomena 
such as the crossover from BCS superfluid[3] to Bose- 
Einstein condensate (BEG) of moleculcs[4], or FFLO 
states in imbalanced Fermi mixtures[5]. Some of the re- 
cent experiments have focused on the case of large popu- 
lation imbalance which leads to the polaronic behavior for 
attractive [6] and repulsive impurity atoms [7]. The results 
are found to be in a good agreement with corresponding 
theoretical treatments[8, 9]. The latter experiments per- 
formed in three-dimensional (3D) traps are commonly in 
a weakly interacting regime where the mean-field uni- 
form theory is valid. On the contrary, experiments 
done in quasi one-dimensional (ID) tubes[10, 11] or 
microtraps[12, 13] have been able to demonstrate un- 
usual phenomena of strongly interacting ID gases, like 
the Tonks- Girardeau (TG) limit [14] where the local cor- 
relation function is totally suppressed 32 = 0[15]. Re- 
cently many theoretical works discuss Fermi-mixtures in 
lower dimensions [16]. In advanced experiments with 
microtraps [12, 13] a few-body ID Fermi-mixture is pre- 
pared with each component consisting of a specific num- 
ber of fermions. Currently [17], the ground-state energy 
of an ideal Fermi gas with a specific number of atoms 
interacting with an impurity was measured with a high 
precision for different interaction strengths. 

In this Letter we provide a theoretical treatment of ID 
trapped ideal Fermi-gas system interacting via a contact 
potential of arbitrary strength with an impurity. We set 
off from the corresponding uniform case for which an ana- 
lytical expression for the ground-state energy was derived 
by McGuire [18], and we analyze the energy contributions 
and the pair correlation function. We demonstrate that 
a modification of this solution offers a precise expression 



to the ground-state energy in the trap, comparing with 
exact numerical results and perturbative methods. Ad- 
ditionally, we perform an analysis of the energy contribu- 
tions and calculate the frequency of the breathing mode, 
which is another important observable for the actual ex- 
periments. Our treatment is directly applicable to the 
ongoing experiments [17] and provides a solid theoretical 
ground for understanding the properties of impurities in 
ID trapped fcrmionic systems or similar physical models. 

The system of a ID ideal Fermi gas with an impurity 
(usually an atom of the same species with equal mass m 
in a different hyperfinc state eg. for ^Li \F ~ 1/2, mp = 
±1/2) [12, 13]) is described by a simple and quite general 
model Hamiltonian: 
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where Ximp is the position of a single impurity, Xi with 
i = 1,A'^ are positions of ideal fermions, and is 
impurity-fermion coupling constant. In the experiments 
in quasi-lD traps[10-13] the transversal trapping fre- 
quency u}± is typically one order of magnitude larger than 
the longitudinal one validating the ID treatment (the 
trap becomes then highly anisotropic with the transversal 
oscillator length a± = yJTiJmu^ much smaller than the 
longitudinal one a\\ = ^Jh/mLO\\). The scattering prop- 
erties of the system are modified by the trapping geom- 
etry leading to confinement-induced resonances[19] and 
arc described effectively by a single parameter, the ID s- 

wave scattering length oid ~ — {l ~ ^^^!/^a °^° ) '*J-/'^3D' 
where 030 is the 3D s-wave scattering length (tuned by 
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magnetic field via Feshbach resonances). Due to anti- 
symmetry, s-wave interactions between identical fermions 
are not possible, while p-wave interactions can be ne- 
glected. Therefore the only relevant interaction in the 
system acts between the impurity and identical fermions 
and can be modeled by a zero-range 5-potential with in- 
teraction strength giD = — 2fi,^/maiD[19]. 

Homogeneous system: In the uniform case, wy = 0, 
the dynamics are governed by the competition between 
the kinetic and the interaction energy. The length scales 
of this problem are the size of the box L, the mean in- 
terparticle distance ~ L/N [or the inverse Fermi 
momentum kp^ = (np)^^] and s-wave scattering length 
aiD- In the thermodynamic limit i — > oo the size of the 
box drops out and the only relevant parameter left is the 
dimensionless Lieb-Liniger parameter 7 = — 2/(paiD)- 

The ground-state energy E' of a homogeneous system 
was found analytically by McGuire[18] in a form of an 
interaction shift AE to the energy of an ideal Fermi gas 
E = NEp/'S + AE with 
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is the Fermi energy. This 
IS valid both for repulsive and attractive 
interactions [18]. We will calculate the local pair cor- 
relation function[15] p2 = (*1'(x)*J„p(a;)#i„p(x)*(a;)) 

(where ^'(x) and ^impix) are fermion and impurity 
field operators, respectively) which is proportional to the 
probability of observing simultaneously the impurity and 
a fermion at position x. In a homogeneous system p2 docs 
not depend on the position and its local value can be 
calculated applying the Hcllmann-Fcynman thcorcm[21] 

dE I dgiY^ = (dH/dgiY)) = P2L leading to 
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mp = (*lmp^"np) = 1/i IS the density of the im- 
purity. As shown in Fig. 1, the pair correlation function 
decreases continuously as the interaction strength 7 in- 
creases. It vanishes in the limit of infinitely strong repul- 
sion 7 -foo /92 — >■ 0, i.e., the impurity avoids meet- 
ing other fermions which mimics Pauli exclusion princi- 
ple. In the opposite limit of strong attraction the pair 
correlation function diverges 7 — ^ — cx3 =^ p2 +00, 
i.e., the impurity forms a zero-size deep bound state 
and its position coincides with the position of one of the 
fermions. In the case of vanishing interactions, 7 = 0, the 
fermions and the impurity are uncorrelated and therefore 
P2 = PPo reduces to a simple product of the correspond- 
ing densities. 

The interaction energy AEint is related to the pair 
correlation function (3) as AEint ~ giT>P2L resulting in 

AEint 
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FIG. 1: (Color online) Thin line: impurity-fermion local pair 
correlation function p2 as given by Eq. (3). Thick lines show 
the interaction shift to the energy in units of Fermi energy: 
red solid line, total energy /S.E, Eq. (2); green dash-dotted 
line, interaction energy AEint, Eq. (4); blue short-dashed line, 
kinetic energy AEkin, Eq. (5). 



The same result can be obtained from the Hellmann- 
Feynman theorem as AEint = gijnidH /dgiu)- On the 
other hand, the shift in the kinetic energy AEkin due to 
the interaction between the impurity and the fermions is 
given by the difference between the total energy (2) and 
the interaction energy (4): 
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Figure 1 shows the dependencies of different contribu- 
tions to the energy shift on the interaction parameter 
7. In the case of strong attractive forces, 7 — > —00, 
the large interaction energy leads to a strong bounding 
between the impurity and a fermion. The interaction 
energy provides the main contribution at intermediate 
repulsive coupling 7 « 3 but, approaching the TG limit 
(7 — > +00), it follows the decrease of the probability 
of contact (^2) and the gas becomes ideal with iV + 1 
fermions, i.e. the impurity is fermionized. 

Trapped system: In the presence of a trap, u)^\ > 0, 
the exact energy is not known, but we will show that 
an expression based on the solution of the homogeneous 
system [Eq. (2)] provides a very precise description. 

First we treat the problem within local density approx- 
imation (LDA). The chemical potential p* of a trapped 
system is approximated as a sum /i* = p + ^mui'^^x^ of the 
homogeneous chemical potential p and the external po- 
tential. The density profile of ideal fermions has a Semi- 
co with the size of 



circular shape plda{x) 



/R2 



the cloud R such that p* = mwy i?^/2. The total number 
of fermions in the trap defines the normalization con- 
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ditions and fixes the chemical potential to /x* = iV/zwy. 
The chemical potential /x* can be interpreted in terms 
of the Fermi momentum of a trapped system kp as 



= nrikpY /2rn with, k^'p = ^2N/a\\. It is interesting to 
note that the energy of N ideal Fermi particles 
obtained within LDA, coincides with the exact result. We 
will look for the energy shift AE = E — Eq to the energy 
of a non-interacting impurity Eq = 1tw\\{N'^ + l)/2. 

A first analytic expression for A£' can be ob- 
tained perturbatively in the case of weak interactions, 
taking for the impurity a Gaussian (non-interacting) 
profile p™p(x™p) = ^7i^exp(-a;f„p/aj^|) and us- 
ing the LDA. The energy shift can be calculated as 

AE = JJ PLDA{x)pimp{Ximp)gm6{x " Ximp)dxdXimp = 

J pLDAix) pimpix)giudx which results in the expression 
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where the modified Bessel functions of first kind Iq, Ii can 
be expanded for large N resulting in a \/N dependence. 

An explicit expression for AE valid for arbitrary 
interaction strength can be obtained by generalizing 
McGuire's expression for the energy (2) of a homoge- 
neous system to the trapped case according to the map- 
ping kp k^p, that is the Fermi momentum of a homo- 
geneous system is substituted by the Fermi momentum 
of a trapped system. Accordingly, we introduce a char- 
acteristic parameter 7* of a trapped system in analogy 
with the Lieb-Liniger parameter 7 = —2'K/{kFCiiT))- In 
a trap the Fermi momentum changes to kp = 



and therefore we define 74 as "ft = 
resulting energy in terms of 74 is 
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We test this expression by confronting it with nu- 
merically exact results of diffusion Monte Carlo (DMC) 
method[22]. This method is based on solving the 
Schrodinger equation in imaginary time and asymptot- 
ically gives the exact ground-state energy. The conver- 
gence is greatly increased by using the importance sam- 
pling according to the guiding wave ipixi, a; at, Ximp) = 



exp(- 



-ax\ 



%rap. 



nN I 



N 

j<k 



Ximp ~ oidIj where parameters a and /3 are optimized 
by minimizing the variational energy. 

The accuracy of the different approaches is analyzed in 
Fig. 2 where dependence of the energy of a trapped sys- 
tem on N for different interaction strength is confronted 
with the exact results of DMC calculations. In the case 
of weak interactions, aiD/a|| = — 10 [Fig. 2a], the modi- 
fied McGuire energy (8) becomes exact for iV ;:g> 1, while 
for few fcrmions a certain difference is present as also 
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FIG. 2: (Color online) Ground-state energy of A'' fermions 
and a single impurity in a trap in oscillator units, (a) weak 
tsiD/a|| = —10 and (b) strong a\Y>la\\ = —1 interaction. Cir- 
cles: DMC points; cross: exact two-particle (one fermion, 
one impurity) energy from Ref. [20]; red solid line: modified 
McGuire expression (8); blue dashed line: perturbative ex- 
pression with as Gaussian impurity -I- LDA , full expression, 
Eq. (6); green short-dashed line: the same, ypN expression, 
Eq. (7). Inset in (b): DMC energy shift versus energy shift 
in modified McGuire expression (8). 



can be seen from comparison with the exact result of 
(8) Busch et al. for N = 1[20]. The perturbative calculation 
agrees well in its full form, Eq. (6), and slightly less accu- 
rate in \/N expansion, Eq. (7). For stronger interaction 
strength, aiD/a\\ = — 1 [Fig. 2b], the modified McGuire 
energy works well for large number of fermions TV. In 
order to quantify the convergence of this expression on 
N we study the correction in the inset of Fig. 2b. 

Once we have demonstrated an excellent accuracy of 
the modified McGuire's formula, we use it to calculate the 
contributions to the total energy. The potential energy of 
the harmonic confinement is related to the mean square 
system size which is calculated using Hcllmann-Feynman 
im(x^): 
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In the limit of vanishing interactions 74 — >■ = 
{N"^ + l)/2 the size of the system is that of a non- 
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FIG. 3: Interaction shift to the energy of a trapped gas with 
an impurity. Black solid line - total energy, Eq. (8), blue short 
dashed line - potential energy of the harmonic confinement, 
Eq. (9), red dash-dotted line - interaction energy, Eq. (10), 
green dashed line - kinetic energy, Eq. (11). 



interacting impurity and ideal fermions. The infinitely 
repulsive (TG) case 7* -> +00 ^ {x'^)/a'^\\ ^ [N + lf/2 
is equivalent to having + 1 ideal fermions. The size in 
the limit of an infinitely strong attraction, 7^ — >■ — cx) 
(x^)/a^l = [N — 1)^/2, is that of — 1 fermions since 
one fermion is absorbed in the bound state with the im- 
purity, which as well can be seen from the separation of 
the energy (8) into binding and ideal Fermi gas energies: 
E = -h^/mal^ + (iV - l)2fiw/2 + 0(7^^)- 

The interaction and kinetic energy can be obtained in 
a similar way as in the homogeneous system: 
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In Fig. 3 we compare the different contributions to the 
energy shift in a trapped system. The interaction en- 
ergy vanishes for 74 = or 74 -I- 00, as the system is 
equivalent to that of ideal fermions, while for 7^ — > — 00 it 
becomes essentially that of a dimer atom-impurity bound 
state. The potential energy contribution to the shift be- 
comes substantial when approaching 74 -> ±00. 

Another experimentally important quantity is the fre- 
quency of the breathing mode w. It is related to the 
response of the system size on changing the frequency of 
the confinement uj^ = -2{x'^)/{d{x^)/duj^^)[2S\: 
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FIG. 4: Square of the frequency of the breathing mode as a 
function of 74 as given by Eq. (12). 



As can be seen from Fig. 4, uj changes significantly, only 
for small number of atoms. In the limit of infinite re- 
pulsion the impurity behaves as an additional fermion, 
leading to the ideal Fermi gas result w^/w^ = 4. Making 
the repulsion finite, the interaction "softens" and the fre- 
quency goes down. In three limiting cases (74 — 0; ±00) 
the result is the same as for ideal fermions making the 
overall dependence non-monotonic. For weak interac- 
tions, 74 — > 0, Gaussian ansatz result (7) is applica- 
ble leading to a linear dependence with 74: w^/o;^ ~ 
4 — Ajt / {'^''^ N) . Also the expansion for strong interac- 
tions 74 ±00 wVwj^ = 4 - 16N/[{N ± 1)^74 - 47V] 
decays as 1/N for large system sizes. 

In conclusion, we provided a theoretical treatment 
for an ideal Fermi gas in a presence of a single impu- 
rity interacting with arbitrarily strong (5-potential. In a 
homogeneous system the total energy was obtained by 
McGuire[18]. Here we analyze potential and kinetic en- 
ergy, as well as experimentally relevant impurity-fermion 
local correlation function[15]. In a trapped system we 
generalize McGuirc's formula and demonstrate in a di- 
rect comparison with DMC result that the obtained ex- 
pression is accurate. This expression in used to calculate 
kinetic, interaction and potential energies. Finally, the 
frequency of the breathing mode is reported for different 
number of fermions. 
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